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Abstract 



We examine various aspects of supersymmetric particle production at lin- 
ear e+e" colliders operating at a centre of mass energy of y/s = 500 GeV, and 
integrated luminosity of / Cdt = 20 — 50 fb^^. Working within the framework 
of the minimal supergravity model with gauge coupling unification and radia- 
tive electroweak symmetry breaking (SUGRA), we study various signatures 
for detection of sparticles, taking into account their cascade decays, and map 
out the regions of parameter space where these are observable. We also exam- 
ine strategies to isolate different SUSY processes from another. In addition, 
we perform four detailed SUGRA case studies and examine the detectability of 
sparticles when several SUSY processes are simultaneously occuring. We show 
that precision mass measurements of neutralinos, sneutrinos and top-squarks 
are possible, in addition to previously studied precision mass measurements 
of sleptons and charginos. 
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I. INTRODUCTION 



The realization that weak scale supersymmetry (SUSY) can stabilize the symmetry 
breaking sector of the Standard Model (SM) has made the search for supersymmetric parti- 
cles one of the standard items for experiments at high energy colliders. Experiments at 
the CERN LEP2 e"^e~ collider will soon probe the existence of charginos, sleptons, squarks 
and even the lightest of the Higgs scalars of supersymmetric models if they are lighter 
than about ~ 80 — 95 GeV . Direct searches for gluinos and squarks lighter than about 
300 GeV will be carried out at the Fermilab Tevatron Main Injector pp collider, which should 
commence operation towards the end of the century. Assuming the unification of gaugino 
masses, experiments at the Main Injector may have their greatest reach for supersymmetry 
(SUSY) via the clean trilepton channel from W1Z2 3i production, and ought to be able to 
probe chargino masses ranging from ~ 50 — 200 GeV (corresponding to gluinos as heavy as 
500-600 GeV), depending on the values of model parameters The supersymmetry reach 
of Tevatron experiments is sensitive to the assumption of i?-parity conservation and may be 
significantly larger (smaller) if the lightest supersymmetric particle (LSP) decays via e or 
H (baryon number) violating interactions 0. It thus appears that while LEP2 or the Main 
Injector could well find a SUSY signal, a decisive exploration of the existence of weak scale 
supersymmetry would require a direct investigation of the TeV scale. This is possible, for 
instance, at the CERN Large Hadron Collider (LHC) which can probe gluino and squark 
masses up to ~ 1300 — 2000 GeV with just 10 fb~^ of integrated luminosity We will 

see later that electron-positron colliders operating at a centre of mass energy of ~ 1.5 TeV 
would have a similar reach for SUSY discovery as the LHC: a SUSY signal may, of course, 
be discovered at a collider operating (much sooner) with a lower value of ^/s. 

If signals for New Physics are first found, either at the LHC or at the first phase of an 
e~^e~ collider operating at 300-500 GeV, the immediate task at hand would be to establish 
their origin. The strategy for this would be quite different at the two facilities. Focussing 
on supersymmetry as the origin of these signals, at the LHC, we would expect signals in 
various event topologies (n-jets plus m leptons plus with and without b tags) from several 
sparticle production processes. It is just this complex plethora of signals that would point 
to their SUSY origin. In contrast, it is unlikely that (except for anticipated degeneracies, 
e.g. for various flavours of sleptons) several SUSY reaction thresholds will be crossed at the 
relatively low centre of mass energy of the first phase of the Linear Collider, so that signals 
may be present from just a single source, and in just a few (perhaps, only one) channels. 
This would, of course, greatly facilitate their interpretation; the price vis-a-vis the LHC is 
that only a few sparticles might be kinematically accessible. Having found indications for 
supersymmetry, the next step would be to sort out the super-particle masses and quantum 
numbers. The sparticle mass spectrum could provide clues about the mechanism of 
supersymmetry breaking, and perhaps also provide a window to the symmetries of GUT or 
Planck scale physics. The detailed measurement of sparticle properties may well be very 
difficult at a hadron collider, owing to, among other things, large backgrounds, indefinite 
subprocess collision energy, additional QCD radiation and the fact that, in general, several 
SUSY reactions may simultaneously contribute to any one signal. 

Many groups have been exploring the physics capabilities of a new, high energy linear 
e~^e~ collider P-|lO[] where the cleanliness of the experimental environment makes it relatively 
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easy to detect signals from the production and subsequent decays of new, heavy particles if 
sufficient luminosity is available. We will refer to this machine as the Next Linear Collider, or 
NLC In its first stage the NLC would operate at a centre of mass energy of ^/s ~ 300 — 500 
GeV, and accumulate about 20-50 fb~^ of integrated luminosity over the first several years 
of operation. It is envisioned that the energy of the machine would be upgraded in stages 
to the TeV region. In addition, there exists the possibility of longitudinal electron beam 
polarization, perhaps reaching magnitudes of ~ 90 — 95% left- or right- polarization. 

The SUSY discovery potential of experiments at such a facility has been the subject of 
many previous studies. In Ref. |T^, it has been shown that using relatively simple cuts. 



WiWi can be detected above backgrounds over almost the entire kinematically accessible 
regions of parameter space by searching for mixed leptonic/hadronic or purely hadronic 
decays of the chargino pair. In Ref. it has been shown that sleptons can usually be 



detected above background for slepton masses up to ~ 90% of the beam energy. Ref. ||T3[ 



shows that with the availability of beam polarization, it should be possible to measure squark 
masses to about ~ 5 GeV, even taking their cascade decays into account: in particular, it is 
possible to obtain the difference between the masses of ql and qn. Finally, in a pioneering 



paper, two very detailed case studies have been performed by Tsukamoto et. al. |jT^ , showing 
that an e+e~ linear collider operating at y/s ~ 300 — 500 GeV can make a variety of precision 
measurements (which test the assumptions underlying the supergravity GUT framework) 
of sparticle masses, spins and coupling parameters. This study makes innovative use of the 
capability for polarization of the electron beam. It is, however, assumed that the beam 
energy is adjustable so that it is sufficient to focus on signals from a single SUSY reaction, 
taken to be chargino or slepton pair production — this simplifies the analysis in that the 
parent sparticle decays directly to the LSP; i.e. there are no complicated cascades to be 
untangled, and further, there is no need to sort out various SUSY processes from one another. 
The purpose of this paper is to expand upon these studies in order to attain two broad goals: 
1.) to delineate the reach of such a machine for supersymmetry and compare it to the reach 
of the LHC, and 2.) to ascertain to what extent precision measurements of sparticle masses 
are possible, even when several SUSY processes are occuring simultaneously, and cascade 
decays of sparticles are operative. 

We work within the framework of the minimal supergravity GUT model (SUGRA) as 



defined in Ref. []T5[ . The only role of "supergravity" is to provide a rationale for the universal 
boundary conditions at an ultra-high scale Mx which we take to be Mqut- This model 
is then completely specified by four SUSY parameters (in addition to SM masses and 
couplings). A hybrid set consisting of the common mass mo {mi 12) for all scalars (gauginos), 
a common SUSY-breaking trilinear coupling Aq all specified at the scale Mx together with 
tan (3 proves to be a convenient choice. These parameters fix the masses and couplings of 
all sparticles. In particular, and the magnitude (but not the sign) of // are fixed. 

The SUGRA framework (and also a SUGRA-inspired MSSM framework without radia- 
tive elecroweak symmetry breaking [|l^) has been incorporated into the event generator 



program ISA JET 7.16 [|1^]. All lowest order 2 — 2 sparticle and Higgs boson production 



mechanisms have been incorporated into ISAJET. These include the following processes ||T8[ 
(neglecting bars over anti-particles): 



e^e qlQl, gnQR, 
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e^e eLCR, 
e+e" ue^e, 

e+e" ^ WiWi, W2W2, W1W2, 

e+e' ^ Z^Zj, = 1- 4), 

e~^e — > ZHg, ZHh, HpHg, HpHh, H^H . 

In the above, £ = e, /z or r. All squarks (and also all sleptons other than stau) are taken 
to be L or i? eigenstates except the stops, for which titi, tit2 and ^2^2 (here, ti,2 being the 
lighter /heavier of the top squark mass eigenstates) production is included. Given a point 
in SUGRA space, and a collider energy, ISAJET generates all allowed production processes 
according to their relative cross sections. The produced sparticles or Higgs bosons are then 
decayed into all kinematically accessible channels, with branching fractions calculated within 
ISAJET. The sparticle decay cascade terminates with the (stable) lightest SUSY particle 
(LSP), taken to be the lightest neutralino {Zi). Final state QCD radiation is included, 
as well as particle hadronization. ISAJET currently neglects spin correlations, sparticle 
decay matrix elements, and also, initial state photon radiation. In the above reactions, spin 
correlation effects are only important for chargino and neutralino pair production, while 
decay matrix elements are only important for three-body sparticle decays. 

To facilitate investigation of polarized beam effects on signal and background cross sec- 



tions 1 19], we have recently incorporated polarized beam effects into the ISAJET e e cross 
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sections. The degree of longitudinal beam polarization has been parametrized as 

Pde') = fi- /i?, where 

, riL 1 + Pl , 

JL = z — — = — — and 

Ir 

In the above, n^^R are the number of left (right) polarized electrons in the beam, and 
fL,R is the corresponding fraction. Thus, a 90% right polarized beam would correspond to 
PL{e~) = —0.8, and a completely unpolarized beam corresponds to PL{e~) = 0. 

In Sec. II of this paper, we make a rough scan of parameter space to delineate the regions 
where a signal for supersymmetry may be observed above background. We also delineate 
those regions where signals from different production processes can be sorted out from one 
other. A complete visual display of the 4+1 dimensional parameter space is of course not 
possible. Instead, we present results in the mo vs. mi/2 parameter plane, for a low and 
intermediate value of tan/3 (tan/3 = 2 and 10), and for both fi < and /i > 0. Variation of 
the parameter Aq mainly leads to changes in third generation squark and slepton masses, 
especially for the top squarks ti and ^2- Hence, we do not address in detail consequences 
of Aq variation, although in Sec. 3 we perform a case study where the ti becomes light 
enough to be accessible. A recent analysis of production and decay rates of third generation 
sparticles at linear e^e~ colliders can be found in Ref. pO|. In most of the work presented 
here, we set Aq = 0. 

Our objective in Sec. Ill can be viewed as a follow-up to reference [|1^]. Here, we perform 
a number of case studies in some detail to explore the capability of an NLC to perform 
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precision measurements of sparticle masses and parameters. We allow for the simultaneous 
production of all sparticles with cross sections and decay patterns as given by the model, and 
describe our attempts to extract the masses of various sparticles. We examine the following 
four cases, which have also been the subject of recent studies on supersymmetry at various 



workshops on linear coUiders ||2T|]. For each case, we list (mo, mi/2, Aq, tan/5 and sign{fi)), 



with the mass parameters in GeV. We have taken nit = 180 GeV throughout this paper. 



(400, 


200, 


0, 2, -1) 


(dominantly chargino production). 


(100, 


300, 


0, 2, -1) 


(dominantly slepton production). 


(200, 


100, 


0, 2, -1) 


(mixed chargino/slepton/sneutrino production). 


(300, 


150, 


-600, 2, 


+ 1) (includes titi production). 



For each case, we use ISA JET 7.16 to simultaneously produce all allowed SUSY particles as 
well as SM backgrounds. 

We conclude in Sec. IV with a summary of our results and some comparisons between the 
NLC and the CERN LHC pp collider. Finally, in an appendix, we collect various expressions 
for lowest order production cross sections of SM and SUSY particles via polarized beams. 



II. REACH OF THE NLC IN SUGRA PARAMETER SPACE 

A. Kinematic Reach 

To gain some orientation for our study of the SUSY reach of the NLC in the parameter 
space of the minimal SUGRA model, we begin by examining where in parameter space the 
various sparticle production mechanisms are kinematically accessible. The lightest SUSY 
Higgs boson Hi is very special, since in the minimal model, rriH^ ^ 120 — 130 GeV p2 



and so, is accessible either via e^e~ ZH^ or via e^e~ HpH^ processes for all values of 
parameters. Furthermore, if is large as is often the case in the SUGRA framework, its 
couplings to SM fermions and gauge bosons are expected to be very similar to those of a SM 
Higgs boson. Recently, detailed studies of the detectability of SM Higgs bosons as well as of 



the Higgs bosons of the minimal SUSY model at the NLC have been performed by Janot P3 
Even for the most difficult case, where rriHi ~ Mz, a 5a signal should be attainable within 
a month of running at or near design luminosity at a 500 GeV linear collider, assuming a b 
tagging efficiency of 50% with a rejection of 97.5% (99.9%) against cc (light quark pairs) is 
achieved. Detection of Hi should be possible even if it decays invisibly via Hi ZiZi pairs. 
Thus, even if the lightest SUSY Higgs boson Hi eludes detection at LEP2, the Tevatron and 
the LHC, it should certainly be discovered at NLC. Hence, the NLC will be able to exclude 



the minimal SUGRA model if no Hi signal is seen |2J] . However, if an Hi signal is the only 



new signal seen, it may well be difficult to distinguish whether one has seen a SUSY 



or a SM Higgs boson: thus, detection of Hi alone will likely not be definitive evidence for 
supersymmetry. In the remainder of this paper, we assume the lightest SUSY Higgs boson 
Hi will be detected at least at the NLC, and so focus our attention on the detectability of 
super-partners. 
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In Fig. 1, we show the regions of the mo vs. mi/2 parameter plane where various 2^2 
SUSY particle processes are kinematically accessible to an e~^e~ linear collider operating at 
a/s = 500 GeV, for tan/3 = 2, Aq = and a) /i < and b) fx > 0. In Fig. 2, we show the 
corresponding results for tan/5 = 10. The regions labelled by TH are excluded by various 
theoretical constraints [|T^, while regions labelled EX are excluded by experimental searches 
for SUSY at the LEP p6| and Fermilab Tevatron colliders as described in Ref. [|l5l (we do 
not include constraints from the LEP1.5 run). The regions below the contours labelled by a 
sparticle pair is where the corresponding reaction is kinematically accessible. The neutralino 
pair contours are an exception. Since the neutralino production cross section can be strongly 
suppressed by mixing angle factors, we conservatively show the regions where the production 
cross section a{ZiZj) > 10 fb (i = 1 — 4, j = 2 — 4). We see from Fig. la and b that the 
outermost boundary of the kinematic reach for SUSY is comprised of three contours: the 
GReR contour for low mo, the WiW^i contour for large mo, and a small intermediate region 
where the Z1Z2 reaction might be accessible. The situation is similar in Fig. 2, except that 
an additional sliver of parameter space may be accessible by searching for fifi pairs (because 
of tl — Tfi mixing induced by the tau Yukawa coupling which increases with tan (3, the lighter 
of the two staus is lighter than cr) as well; in practice, the detection efficiency is larger 
for identifying selectron and smuon signals so the r channel is unlikely to be relevant for 
the maximal reach. The smuon contours (as well as stau contours in Fig. 1) essentially 
overlap with the selectron contours and have not been shown. However, the detection and 
measurement of their properties is very important for testing slepton universality |l^ and 
the flavour structure of the slepton sector: this could shed light [^] on the nature of physics 
at the GUT scale. Interesting information about the composition (gaugino versus Higgsino) 
of Zi 1^ can also be obtained from a detailed study of the stau signal. For parameter space 
points outside these regions, SUSY particles will be accessible at NLC500 only via higher 
order processes p9| , while inside these regions, at least one and often many SUSY particles 
might be produced with significant rates. 

In the low mo region, we see that various slepton pair reactions- d-Rd-R, eRCi, dhd-L and 
V(yr- can be accessible. The sneutrinos may or may not have visible decay products, de- 
pending on how massive they are relative to the charginos and neutralinos. Potentially, 
the greatest reach for is via the tRtL production which occurs via t— channel neutralino 
exchange; the production cross section is, however, significantly smaller ||18[ than cr{eReR) 
or (7(6^6^) (if both selectrons have the same mass). Smuons and staus can only be produced 
via s-channel processes. The regions with mo ^ 250 GeV are also the regions most favored 
by cosmological neutralino relic density constraints [|30|, which require 100 < m^ < 250 GeV 
to obtain a dark matter relic density in accord with infiationary cosmological models with 
1:2 mixed hot to cold dark matter. 

In the large mo region, sfermions are too heavy so that this part of parameter space can 
best be probed by searching for chargino pairs, and possibly Z1Z2 or Z2Z2 pairs. Since Zi 2 
is gaugino-like for a wide range of SUGRA parameters, these neutralino pair cross sections 
depend significantly on t or u channel or cr exchange, so the neutralino pair production 
rate typically decreases with increasing mo. The large mo region is difficult to accommodote 
cosmologically since the annihilation rate of the (stable) LSPs via sfermion exchanges is 
suppressed; this implies too short an age for our universe unless ZiZi annihilation can 
efficiently proceed via s-channel Z or resonances 
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In the regions with smaller mo and mi/2, the higher mass chargino and neutralino states 
W2, Z3 and Z4 may become accessible to NLC500 searches. For still lower mo and mi/2, iiii 
production and ultimately qq production can become accessible. However, for NLC500, rriq 
must be less than 250 GeV to be at all accessible, so that these strongly interacting states 
would have been seen much earlier at Tevatron collider experiments. Finally, we see that 
in the lowest mo and mi/2 regions, heavy Higgs bosons such as may become accessible 
to NLC searches. 

By comparing Fig. la and b with Fig. 2a and b, we see that although the sparticle 
accessibility contours can change somewhat with variations in tan/? or sign{n), the overall 
qualitative trends are the same: NLC500 is most likely to access the non-colored, charged 
SUSY particle states and sneutrinos, with little hope of seeing squarks or gluinos, unless 
they have already been discovered at Tevatron experiments. The exception is the possibility 
of seeing third generation ti or bi squarks, which can have significantly lighter masses than 
the other squarks, which should be nearly mass degenerate in the minimal SUGRA model. 
Variation of the Aq parameter mainly affects third generation ti, bi and fj masses (i = 1 or 
2), so that only their reach contours change appreciably with Aq. 

B. Results from Event Simulation 

We now turn to the evaluation of the reach of the NLC by comparing signal against SM 
backgrounds using explicit event generation. Since beam polarization has been shown to be 
a useful tool, as a first step, we show in Fig. 3 various lowest order SM background cross 
sections, as a function of electron polarization P/;,(e~) for an unpolarized positron beam. The 
e~^e~ —>■ e+e~ contribution includes only s-channel contributions. For unpolarized beams 
(PL(e~) = 0), WW production is the dominant SM process. By tuning PL{e~) to ~ —0.9 
(95% right polarized beam), the magnitude of the WW cross section can be significantly 
reduced relative to other SM backgrounds, which show only a mild dependence on beam 
polarization. Since WW production is a major background for many new physics processes. 
Fig. 3 suggests that the use of (dominantly) right-handed electron beams would yield a better 
signal to background ratio, except for those signals whose cross sections become small when 
Pi[e~) ~ — 1. We have not shown backgrounds from 2 — 3 and 2 — 4 SM processes- these 
can be reduced by using suitable cuts 

Our next step is to generate explicit events for signal and background. We focus on 
optimizing cuts for eRCR and WiWi production, which should give the largest reach into pa- 
rameter space. We use the ISAJET toy detector ISAPLT with the following characteristics. 
We simulate calorimetry covering —4 < 77 < 4 with cell size Arj x A0 = 0.05 x 0.05. En- 
ergy resolution for electrons, hadrons and muons is taken to be A_E = \J .Q22hE + (.01-E)^, 

Ai? = ^J.IGE + {.03Ey and Apx = 5 x lO'^^p"^ respectively. Jets are found using fixed cones 
of size R = a/A?7^ + A0^ = 0.6 using the ISAJET routine GET JET (modified for clustering 
on energy rather than transverse energy). Clusters with E > 5 GeV and |r7(jet)| < 2.5 are 
labeled as jets. Muons and electrons are classified as isolated if they have E > 5 GeV, 
\rj{i)\ < 2.5, and the visible activity within a cone of R = 0.5 about the lepton direction is 
less than max{^, 1 GeV). Finally, 6-jets are tagged with an efficiency of 50%, while c-jets 
are misidentified as 6's with an efficiency of 3%. Jets with one or three charged prongs are 
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classified as rs for tlie purpose of r-veto (see Sec. HID). 

The signature for IrIr production is a pair of acollinear same flavor/opposite sign leptons 
recoiling against I^. To search for such a signal, we essentially follow the cuts of Ref. 
and require i) 5 GeV < E{1) < 200 GeV, ii) 20 GeV < E{visihle) < 400 GeV, in) 
|m(£+r) - Mz\ > 10 GeV, iv) \ cos^(£±)| < 0.9, v) -QecosOe < 0.75, vi) Oacop. > 30°, viz) 
I^rp > 25 GeV and viii) veto events with any jet activity, where the polar angle is measured 
from the electron beam, and Q is the charge of the lepton. Cut in) eliminates backgrounds 
from e~^e~ — ^ ZZ^vvZ and e^e~Z production, while cuts iv) and v) greatly reduce the 
backgrounds from WW and evW production (we neglect the latter). For unpolarized beams, 
the resulting background level was 17 fb, while for Pi^[e~) = —0.9, the background was 2.4 
fb. Thus, for the polarized case, a 5a signal for 20 fb~^ of integrated luminosity requires a 
signal rate larger than 1.73 fb. 

To search for chargino pairs, one may search either for 4-jet events from both Wi hadronic 
decays, or l-£+2-jet events from mixed hadronic/leptonic chargino decays. We found that 
either signature gives a similar reach; we focus on the mixed hadronic/leptonic signature 



since ultimately it is more useful for chargino mass measurements. Following Ref. |T^, we 
require events with one lepton plus 2 jets, and i) # of charged tracks > 5, li) 20 GeV 
< E{visible) < 400 GeV, ni) if E{jj) > 200 GeV, then m{jj) < 68 GeV, %v) > 25 GeV, 
v) \m{iiy) - Mw\ > 10 GeV for a H^-pair hypothesis, vi) |cos6'(j)| < 0.9, |cos6'(£)| < 0.9, 
-QiCos6c < 0.75 and QiCos6{jj) < 0.75 vii) OacopXWW) > 30° for a W^-pair hypothesis. 
Although the dominant background from WW production is smallest for Pii^e^) close to -1, 
the signal cross section also drops rapidly since the chargino is frequently an SU{2) gaugino, 
and so, couples only to the doublet electron. Hence, the use of left-handed electron beams is 
required. For Pi(e~) = +0.9, the resultant background level was 155 fb, so that a 5cr signal 
for 20 fb~^ of integrated luminosity requires a signal rate larger than 14 fb. 

In Fig. 4 and Fig. 5, we show the 5cr reach of NLC500 for minimal SUGRA via 
and l-£ + 2-jets searches for tan (3 = 2 and tan (3 = 10, respectively. Here, we assume an 
integrated luminosity of 20 fb~^ and compare the reach we obtain with the kinematic reach 
contours of Figs. 1 and 2, shown as dashed contours for the and li + 2j signals. In 
Fig. 4a alone, we compare the reach for a polarized e~ beam with the unpolarized case. 
The dotted curves correspond to the NLC500 reach using the above cuts with unpolarized 
beams, while the solid curves correspond to the reach for a 95% polarized electron beam 
with dominantly left (right) handed polarization for the chargino {£r) search. Notice that 
for mo ^ 250 GeV, the signal (with unpolarized beams) from charginos is observable 
in between the two dotted curves; below the lowest dotted curve the chargino is rather 
light, and our cuts are not optimised for their selection. By comparing the dotted and 
solid curves, we see that there is only a marginal gain in reach using polarized beams with 
PL{e~) = —0.9 for the slepton signal (to reduce WW backgrounds), and PL{e~) = +0.9 for 
the chargino signal (to gain the largest signal cross-section). For this reason, we have chosen 
not to show the polarization dependence in the other frames in Fig. 4 and Fig. 5. The real 
power of polarization is for precision measurements of masses and couplings [H,^. We 



note that the selectron search contours fill most of the region of slepton accessibility, except 
for a small region around (mo,mi/2) = (100,500) GeV, where m^^ ~ ""^Zi' 
leptons have very little visible energy. The £+ 2-jets signal from chargino pair production 
can likewise be seen almost up to the kinematic limit over much of the parameter space; 
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the exception is around {mo, 1711/2) = (250,275) GeV, where the i>L becomes hght enough to 
cause a significant drop in the chargino pair production total cross section. For even lower 
values of mg, the region of detectable charginos falls off due to rising (diminishing) chargino 
leptonic (hadronic) branching fractions, as sleptons become very light, in which case a signal 
may be observable in the acollinear efi + I^j, channel. 

Some additional reach may be gained by looking for Z1Z2 production around 
(mo,mi/2) = (250,320) GeV. For the tan/? = 2 cases, B{Z2 ZiHe) is almost 100%, 
so the Z1Z2 signature is bb + If. The physics background consists of ZZ and ZHi produc- 
tion, which occurs at a much larger rate. For the tan/3 = 10 cases, Z2 Z^Z is comparable 
to Z2 ZiHi, so one can also look for Z+If events, where Z i'^i^ . The Z+If signals also 
have large backgrounds from WW production. To determine the viability of the Z1Z2 sig- 
nal in this region, we examined the point {mo,mi/2, Aq, ta.n (3, sgn{fi) = (240, 340, 0, 2, +1) 
in parameter space. In our simulation we used Pl{g~) = +0.9, and assumed 20 fb~^ of 
integrated luminosity. To isolate the Z1Z2 signal, we require events with 2 tagged 6-jets, 

> 25 GeV and 30° < Acpi^i < 150°. We also require that the missing mass > 340 GeV 
since for this case, the two main backgrounds ideally have = Mz, up to considerable 
smearing corrections; in contrast, the SUSY signal requires > 2m^^. For the particular 

point we considered for our simulation, we found a total Z1Z2 cross section of 21 fb, with a 
signal efficiency of 6%. No SUSY or SM background events were found. We thus conclude 
that a left-polarized cross section of a{ZiZ2) ~ 10 fb is needed to achieve a ~ 10 event signal 
with 20 fb~^ of data. We show in Figs. 4 and 5 the 10 fb contour labelled Z1Z2, which 
represents the rough reach in parameter space for the Z1Z2 signal, assuming that Z2 —>■ ZHi 
is the dominant decay of Z2 and that the detection efficiency varies slowly in this region of 
parameter space. 

III. SPARTICLE MASS MEASUREMENTS AT LINEAR COLLIDERS: FOUR 

CASE STUDIES 

Our task in this Section is to go beyond the detectability of new particles, and to examine 
prospects for the determination of their properties. Towards this end, we perform four 
detailed case studies where we attempt to isolate the different SUSY production processes 
from one another in order to facilitate the interpretation of our results. We focus here on 
precision mass measurements, although certainly a wide range of other measurements such as 
spin quantum numbers, total cross sections, sparticle branching fractions etc. are possible 
||14|| . Such measurements can serve as tests of the underlying framework (the minimal 
SUGRA model, in our case), and perhaps even help to determine some of its fundamental 
parameters. 

Since SUSY particle decays always terminate in the LSP Zi, a direct reconstruction of 
SUSY particle masses via "mass bumps" is not possible. However, the cleanliness of e^e~ 
scattering events, combined with the well-defined initial state, leads to kinematic restrictions 
which depend directly on sparticle masses. For instance, in the reaction e^e~ —^pi+ P2, 
followed by P2 ^ Ps + Pa, the energy of particle p^ is restricted to lie between 

7(^3 - m) <E3< 7(^3 + m), (3.1) 
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where 



(ml + mj - ml)/2m2, pI = \/ 



- ml 
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1 — 1/72 and 



E2 = {s + rrig — m^)/2y^, up to corrections from energy mis- measurements, particle losses, 
bremsstrahlung, etc. We will see below that this formula provides a simple yet clean way 
for the determination of slepton and LSP (or sneutrino and chargino) masses |jl4| and, with 
appropriate analysis, also of the chargino mass when the chargino decays via Wi —>■ ffZi 
(see Sec. HID). 



A. Case 1: Dominantly chargino production 

The first case that we examine corresponds to the SUGRA point 
(mo, mi/2, Aq, tan (3, sgn{fi)) = (400, 200, 0, 2, —1), (where parameters with mass dimensions 
are in GeV) whose locus in the mo vs. mi/2 plane is labelled "1" in Fig. 4a. We show in Fig. 
6 the total cross sections for accessible 2 — >^ 2 SUSY particle reactions as a function of beam 
polarization PL{e~). These may be compared directly to various SM cross sections shown in 
Fig. 3. We also plot the masses of the accessible SUSY particles to help orient the reader. 
For this case, m^^^ = 85 GeV, so it would likely have been already discovered at LEP2. The 
ZHg cross section only varies by 50% over the range of PL{e~). The dominant SUSY process 
is lyiW^i production {rriy^^^ = 175 GeV) and, because the chargino is essentially an SU{2) 
gaugino, its cross section drops rapidly as Piie^) —1. At high energy (i/i >> Mz), we 
may think of just the neutral SU{2) vector boson exchange contributing in the s-channel, 
so that s-channel amplitudes for right-handed electrons are suppressed; since the sneutrino 
exchange amplitude always involves just left-handed electrons, the polarization dependence 
of this cross section is readily understood. The same reasoning explains the behaviour of 
the Z2Z2 cross section. In the limit that Zi is the bino, the t-channel selectron amplitude 
dominates ZiZi (and Z1Z2) production. The polarization dependence of a{ZiZi) is readily 
understood once we recognize that the cross section varies as V^, where Y is the hypercharge 
of the selectron exchanged in the t-channel. Finally, because Z2 has suppressed hypercharge 
gauge couplings, the polarization dependence of a{ZiZ2) follows that of a{Z2Z2). 

In this scenario, Wi WZi with nearly 100% branching ratio, so the WiWi signal should 
be easily seen above the 5a level of Sec. 2 in either the 4-jet or l£+2-jets mode. However, 
to extract a chargino mass, a clean event sample is needed, and further discrimination of 
signal from SM (mainly WW) background is necessary. We focus here on the l£+2-jets 
signal, for which mass measurements are relatively straightforward. We use unpolarized 

beams and assume / Cdt = 20 fb~^. The missing mass, defined by r?^ = provides 
a powerful discriminator. For WiWi production, is constrained to be > 2m^^ = 172 
GeV, while WW production has no such constraint. We show in Fig. 7a the distribution 
for both signal and background. In this case, a rather clean SUSY signal can be obtained 
by requiring > 240 GeV. The distribution of surviving events is plotted as a function of 
dijet energy Ejj in Fig. 7b. The background level is indicated by the histogram, while the 
signal cross section is shown by the points with error bars. In this case, Ejj ~ E^ from 
the Wi — > WZi decay, so the endpoint structure of a 2-body decay discussed above should 
apply. The tips of the arrows indicate the the theoretically expected endpoints obtained 
using Eq. (3.1). The distribution has significant smearing (particularly at the low end) due 
to our calorimeter simulation and use of the cone algorithm for jet finding (which entails 
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some loss from energy outside the cone), and may well be improved with different jet finding 
schemes, or by plotting all visible energy aside from the detected lepton. Since Tsukamoto 
et. al. have already shown that a fit to the Ejj distribution leads to a mass measurement 
of to ~ 5%, we have not made any attempt to improve our jet algorithm and repeat 
this same analysis here. 

Instead we focus on the other interesting possibility which is to isolate a signal from 
Z2 in order to measure its mass. For our case, B{Z2 ZiHi) = 99.6%, so that Z1Z2 
production almost exclusively results in bb + events. The physics background here is 
mainly due to ZZ and ZHi production. In this case, we use Pi(e~) = +0.9, and assume 50 
fb~^ of integrated luminosity. To isolate the Z1Z2 signal, we require events with two tagged 
6-jets, I^j, > 25 GeV and 30° < A^^g < 150°. At this point, we can proceed with a plot 
of n^, which is shown in Fig. 8a. For this case, the two main backgrounds have = Mz, 
up to considerable smearing corrections, while the SUSY signal requires > 2m^^ = 172 
GeV again. In fact, the kinematic endpoints for in the SUSY case are easily calculable. 



since njb = \J if' — ^'^ = \J s — ^^Eh^ + "^Hi (with bounded as in Eq. (3.1)), and are 
indicated on the plot. A clean separation between signal and background can be obtained 
by requiring vfji > 300 GeV. If the Hi mass is well measured from LEP2 or NLC, then these 
endpoints can be used to measure m^^ and m^^. The dijet invariant mass from Hi — >■ bb may 
be checked for consistency with a Higgs boson mass hypothesis. We show its distribution 
in Fig. 8b after the > 300 GeV cut. It peaks somewhat below ttihi due to imperfect jet 
reconstruction, energy mismeasurement and energy loss due to neutrinos. Finally, Eh^ = E^i 
may also be used for a determination of m^^ and fTf^^- We show the dijet energy distribution 
in Fig. 8c along with the expected background (solid histogram). The tips of the arrows 
denote the theoretically expected end points. Once again, we see that there is considerable 
smearing at the lower end. 

The missing mass distribution in Fig. 8 a appears best suited for a mass measurement 
because missing energy from mismeasurement or calorimeter losses cancels out to some 
degree in constructing rrji. We perform a fit to the rr^ distributions, which depends on 



m 



and m^^ (assuming is known from LEP2 or NLC), and plot the resulting Ax^ 



~ Xmin = 2.3 and 4.6 contours (these correspond to 68%, i.e. "la" and 90% CL error 
ellipses) for a 50 fb~^ "data" sample. The result is shown in Fig. 9a, from which we see that 
m^^ = 176.5 ± 4 GeV, and m-^^ = 86.1 ± 3 GeV (Icr error bars), to be compared with the 
input values of m^^ = 175.2 GeV and m^^ = 85.9 GeV. The "data" (points) and the best 
fit (histogram) are shown below in Fig. 9b. 



B. Case 2: Dominantly selectron production 

To contrast with Case 1, we now consider a set of parameters for which selectron pair pro- 
duction dominates, and pair production of charginos is kinematically forbidden at NLC500. 
For Case 2, we choose {mo,mi/2, Aq, tan P, sgn{fi)) = (100,300,0,2,-1). The cross sections 
for accessible processes are shown versus PL{e~) in Fig. 10 together with relevant sparticle 
masses. The polarization dependence of the slepton pair production cross sections is straight- 
forward to understand — in pair production is maximized for a right-handed electron beam 
because it occurs essentially due to hypercharge gauge boson exchange, and Y^Er) = 2Y[(.l) 
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(for enCR production, this effect is further accentuated by the fact that the t-channel ex- 
change amphtude vanishes for Pf^(e~) = 1). The production of clCr pairs, which occurs 
only via t-channel exchange, is independent of beam polarization. Finally, the production of 
left slepton (or sneutrino) pairs shows the opposite dependence as pair production, but 
it should be kept in mind that these have both hypercharge and SU{2) couplings so that 
the analysis for this reaction is not as simple. 

For the parameter choice in Case 2, = 93 GeV and is very near the limit for 

observability at LEP2. We see from Fig. 10 that eRen is the dominant SUSY particle cross 
section over most of the range of Pi^e'). v^vi production also occurs, but since z/£ — > ZxVii 
this process is invisible. e^fiL also occurs at a large rate and can provide an opportunity 
for mass measurement [pieL^i) is an order of magnitude smaller). The cross sections 
for IjiiR {^l^l) production shown in the figure are summed over fi and f, and are smaller 
than the corresponding selectron production cross sections because there is no t-channel 
contribution in their case. Finally, we note that — >■ iZi and are the only 

allowed two body decays of sleptons, so that the analysis of slepton production which is 
signalled by like flavour, opposite sign acollinear lepton pairs, is free from complications 
from cascade decays. 

Observation of a signal in the acollinear e~^e~ + ^ or fi~^fj^~ + channels without an 
accompanying signal in the e'^fi^ channel would suggest a slepton hypothesis. Focussing, 
for the moment on the dimuon channel, we can select out fj,Rp,R production by adjusting the 
beam polarization. Operating with Pi^e') = —0.9 reduces the WW fifi background to 



tiny levels, while enhancing production of right sleptons ||1J]. Turning to prospects for the 
smuon mass measurement, we show in Fig. 11a the E^^ distribution in dimuon events after 
the dilepton cuts of Sec. 2, for an integrated luminosity of 20 fb~^. The background level 
is denoted by the dashed histogram while the arrows denote the theoretical endpoints. In 
Ref. fl^, it was shown a mass measurement of m^^ and can be made to ~ ±1%, so we 
do not repeat the analysis here. 

In Fig. lib, we show the Ee+ distribution for the same luminosity and polarization. 
This distribution consists of two components: the solid histogram is from enCR production, 
while the dashed histogram is from e]^e~^ production. (Note: with PL{e~) = —0.9, e^ej 
is produced at a much higher rate than e^e^.) Since rrig^ = 238 GeV here, in contrast to 
rrien ~ ^^"^ GeV, the endpoints for the two components differ significantly. Now, because 
we already know from Fig. 11a, the upper endpoint will yield mg^ to a similar degree 
of precision (a one-parameter fit may be made). The energy distribution of e~ is shown 
in Fig. 11c, again for both e^eij and eRCi components (notice the different scale from 
Fig. life). The endpoints of 6^6 l are completely contained within the eRCR distribution. A 
clean sample of crCl events can be isolated by requiring Ef.+ > 75 GeV from Fig. lib. Then 
the subsequent E^.- distribution can be plotted in Fig. lid. These endpoints now depend 
on mg^ , rrig^ and m^^ , so that consistency with the previous mass measurements as well as 
with the assumed production mechanisms, can be checked. 

The other opportunity, in this case, is to identify Z1Z2 production and measure fn-^^. 

In this case, Z2 Hl iiZi with a branching fraction of about 5%, so one can search 
again for acollinear dilepton pairs. SM backgrounds mainly come from WW production 
and we have to discriminate the Z2 signal from i£ production processes. We focus here on 
the dimuon signature, so we can ignore background from clCr production. In this case. 
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the fiR and Zi masses are already well measured. We run with PL{e~) = +0.9, for a 
50 fb~^ data sample. To eliminate the WW background, we require > 250 GeV, and 
A(f){fifi) < 90°. In the remaining event sample, the muons from fti decay are very hard, 
so we require E^{fast) > 75 GeV, and plot the E^{slow) distribution, which is shown in 
Fig. 12. A small remaining background from /XL/tL production populates the E^{slow) > 60 
GeV region, while the residual SM background populates the E^{slow) < 20 GeV region. 
The Z1Z2 signal gives a distinct upper endpoint, which depends on Zi, Z2 and fli masses. 
By combining this information with the previous mass measurements, a constraint on m^^ 
can be obtained. 



C. Case 3: Mixed chargino, slepton and sneutrino production 

In order to study the additional complications that arise when charginos and sleptons are 
simultaneously accessible, we are led to consider Case 3 with {mo, mi/2, Aq, tan (3, sgn{fi) = 
(200,100,0,2,-1), for which the spectrum of sparticles and production cross sections are 
shown in Fig. 13. The polarization dependence of the various cross sections is as in Fig. 6 
and Fig. 10 and needs no further discussion. 

Clearly, many superpartners would be produced at NLC500. In this case, rriHe — 69 
GeV, so that the light Higgs boson would presumably have been discovered and studied 
at LEP2. As can be seen from Fig. 13, Wi, Z2, ^l, in, vl and possibly even Z3, Z^ and 
W2 should be accessible. From the cross sections shown in Fig. 13, we see that for right 
polarized beams, CRfiR is the dominant process, while for left polarized beams, in fact VeVe 
is dominant, followed closely by WiWi production. In this case, the sneutrinos and sleptons 
are significantly heavier than the charginos and neutralinos, so that their cascade decays 
need to be incorporated in the analysis. For instance, the u with a mass of 207 GeV, decays 
visibly via u —>■ Z2V (32%), and v — > W\l (61%), with v — >■ Z\v making up the balance. The 
decay pattern of is similar although its direct decay to the LSP occurs about 20% of the 
time. In contrast, Ir which has no SU{2) gauge interactions essentially always decays via 

—>■ iZi. As can be seen from Fig. 13, there are so many SUSY processes taking place at 
\/s = 500 GeV that it is potentially difficult to isolate one from another! The ability to tune 
the beam energy would certainly be desirable in this situation, so that one could sequentially 
study each SUSY process as one passes production threshold. With — 

~ 96 GeV, 

however, at least WiWi, Z1Z2 and Z2Z2 would be occuring simultaneously, even running 
NLC at energies as low as ^/s ~ 300 GeV. Incidently, it is worth remarking that for this 
scenario, SUSY would certainly have been discovered at the Main Injector upgrade of the 
Tevatron, and in several channels, so that it is again reasonable for us to focus our attention 
on precision measurements of sparticle properties. 

With unpolarized beams, one may attempt to measure and m^^ via the WiWi — *• 

iuZiqq'Zi mode. We adopt the chargino cuts of Sec. 2, but find substantial background from 
z>iZ>i production which distorts the Ejj distribution. In this case, running NLC at 400 GeV, 
below ulUl threshold allows a clean distribution of Ejj to be made. A mass measurement 
of '^fy^ and m^^ should be possible with a precision of a few percent [|T^ . 



As a second measurement, we run with PL{e ) = —0.9 in an effort to pick out a eRCR 
signal. We run with the slepton cuts of Sec. 2, for 20 fb~^. The resulting distribution for 
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SUSY signal as well as the SM and other SUSY backgrounds is shown by the histograms 
in Fig. 14. The SUSY background comes mainly from crCl production, and has a slightly 
different endpoint structure. A mg^ and m^^ mass measurement to a few percent should be 



possible from this distribution, as shown in Ref. |]14 . 

The large production cross section for VeVe events for left polarized beams can lead to 
some unique, almost background free signatures. For instance, if one sneutrino decays via 
i>e — >■ eWi — > e+gg'Zi and the other decays via z>e eWi —>■ e+fiuZi, we obtain spectacular 
ee/i + jets + If events. For this topology, we require in addition > 25 GeV and > 40 
GeV. We then plot the distribution, which should have endpoints determined by m,>^ and 
m^^. We fit an appropriate function to the expected distribution, and map out the 
values in the vs. iti^^ plane for a data set of 20 fb~^ at Pi{e^) = +0.9. The minimum 

is shown in Fig. 15a, along with contours of Ax^ = 2.3 and 4.6 from the minimum. 
We obtain a measured value of m^^ = 207.5 ± 2.5 GeV and m^^^ = 96.9 ±1.2 GeV- a 
1% measurement of these masses. The Eg distribution from data, along with the best fit, 
are shown in Fig. 156. Because our fit to the Ee distribution had been done for the same 
"theory" set, we have double checked our procedure by repeating it for a nearby point in 
parameter space. The results are shown in Fig. 15c and 15 d; we see that we obtain similar 
precision as in Fig. 15 a. 

The cascade decays of selectrons and smuons can also lead to unique event signatures 
from eie-L or p,Lp'L production (of course, selectron production dominates). One possibility 
is where each ei —>■ eZ2 — > e + ii + Zi, which can give events with 6 leptons plus If. We 
focus only on this unusual channel, although similar or perhaps even better measurements 
can possibly be performed in other channels. We require events with six leptons, no jets 
and Ifrp > 25 GeV. We further require the two fastest leptons to be of same flavor /opposite 
sign, and then plot Ee of the two hardest leptons. We study a 50 fb~^ sample with Pi^e') = 
+0.9, and obtain a sample with 70 events, all from the signal (this implicitly assumes that 
the left selectrons and smuons are mass degenerate). The energy distribution of the two 
fastest leptons has endpoints depending on m^^ and m-^^. Again, by fitting this distribution, 
contours of Ax^ values are calculated, with the result shown in Fig. 16 a. We see that 
a measurement of m^^ = 221.6 ± 6 GeV and m^^ = 94.7 ± 6.5 GeV is possible. The 
corresponding data and best fit are shown below in Fig. 16b. We note that the relationship 
m| —ml = — cos 2/3 depends only on the SU{2) gauge invariance and so is very robust. 
A good measurement of m^^ and rriu can lead to a model-independent determination of the 

parameter tan f3. Unfortunately for the present case, £l and i) are very close in mass and a 
combination of these mass measurements only implies tan/? > 1.8. A better determination 
of tan j3 may be possible for other values of parameters. 

Finally, for Case 3, we examine Z1Z2 production. We require the dilepton cuts of Sec. 
2, but in addition require 2 opposite sign muons with E^s. < 200 GeV and 4>{ii) < 90°. 
We show the resulting dimuon invariant mass distribution in Fig. 17. The solid histogram 
denotes the "signal" and the dashed one the background. In this case, the signal consists 
of 68% Z1Z2 production, 21% z/z/ production and 10% Z2Z2 production. All these signal 
processes lead to leptonically decaying Z2 plus missing energy. The m^ji'^ n^) plot has a 
relatively sharp upper cutoff at m^^ — m^^ , allowing an independent determination of m^^ 

given the information on from or Wi production processes discussed above. 
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D. Case 4: Includes signals from top squark production 



Up to now, we have only considered cases with Aq = 0. The weak scale A-parameters, 
of course, do not vanish, and are obtained by renormalization group evolution. The A- 
parameters (and hence Aq) mainly affect the phenomenology of third generation sfermions 
(and of gluinos, charginos and neutralinos via modifications of their decay patterns p2| ) 
which can be significantly lighter than their first and second generation siblings. We are thus 
led to consider Case 4 with (mg, mi/2, Aq, tan/3, sgn{fi)) = (300, 150, —600, 2, +1), chosen so 
as to lead to a sparticle mass spectrum which contains a light top squark with m^^ = 180 
GeV which essentially always decays via ti bWi. The SUSY particle masses and total 
cross sections as a function of Piie') are shown in Fig. 18. In this case, rriHt = 102 GeV 
and = 110 GeV, so both are just beyond the reach of LEP2. The novel feature, not 
encountered previously, is the accessibility of iiii production with a cross section of about 
40 fb, independent of PL{e~). A stop with this mass is unlikely to be observable even at the 
Main Injector upgrade of the Tevatron 



Charginos are the most copiously produced sparticles if the electron beam is unpolarized 
or if its polarization is dominantly left-handed. We first consider the prospects for measuring 
Wi and Zi masses by analyzing WiWi production, using^Pi(e~) = +0.9. The difference 
from the chargino mass analysis in Case 1 is that now Wi decays into 3-body qq'Zi and 
iuZi final states (with essentially the same branching fractions as for W decay) so that 
Eq. (3.1) is not applicable. To isolate the WiWi signal from SM background, we use the 
l£+2-jets cuts of Sec. II, and in addition require > 240 GeV. The resulting distribution 
in Ejj is shown in Fig. 19a. The lower end point of this distribution depends on the jet 
algorithm, while the upper end point is given by Eq. (3.1) with = rrijj = 0. The 
mass analysis is thus not as straightforward as in the earlier cases. In this case, we first 
consider the scatter plot of Ejj vs. rrijj shown in Fig. 196. For each value of mjj, Eq. (3.1) 
implies a definite range of Ejj. The resulting envelope of the {rrijj dependent) endpoints of 
the Ejj distribution is shown as the solid contour in the Figure. The chargino signal lies 
almost entirely within this envelope, while some remaining background events populate the 
outer regions. For an integrated luminosity of 50 fb~^ the event sample is large enough to 
artificially force 2-body kinematics onto the Wi decay by requiring that the dijet invariant 
mass rujj lie within some specified bin. In our analysis, we take four bins each of width 4 
GeV, centered at rrijj = 22, 26, 30 and 34 GeV. For each bin, the Ejj distribution follows 
the form for Wi ZiW*, where Mw* is given by one of the previous bin central values. 
The procedure for fitting rn^^ and m^^ for the forced two-body kinematics is the same as 
in previous cases, except that now we fit simultaneously to four different mass bins; the 
resulting contours of Ax^ = ~ Xmin shown in Fig. 20a for a data set of 50 fb~^ 
We find that = 110.6 ± 5 GeV (a 5% measurement), compared to the input value of 
rUy^^^ = 109.8 GeV. Likewise, m^^ = 57.5 ± 2.5 GeV, compared with the input of m^^ = 57 
GeV. These measurements may be improved by trying different bin choices or different 
fitting procedures. The Ejj distribution and best fit to the data for the rrijj bin centered at 
30 GeV are shown in Fig. 20b. In this mass bin, the function 

FiE) = N{1 + exp\{E^,.^ + 21.5 - E)/3.7/a^_J)-i 
x(l + exp[{-Ema. + 24 + E)ll.^laE^jr^ 
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(where all energy parameters are in GeV) provides a fit to the theoretical expectation for 
the Ejj distribution. 

Next, we focus on Z2Z2 events with a Pii^e^) = +0.9 beam. The Z2 here decays via 
the three-body mode, with B{Z2 eeZi) = 4.5%. To obtain a clean event sample, we 
require four isolated leptons with no jets in each event, and require as well I^rp > 25 GeV 
and 20 GeV< Ems. < 400 GeV. We look at only e^e~fi~^fi~ events, and veto any events with 
m{i^i~) = Mz ± 10 GeV. We are left with just 9 signal events for / Cdt = 50 fb~^, with no 
SM background. The distribution of the two like-flavour dilepton masses in each event is 
shown in Fig. 21. This dilepton invariant mass is restricted to lie between and m^^ —m^^, 
so that the highest mass value gives a lower bound on this mass difference (the fact that we 
have a small number of events only means that there is a significant chance that we may 
not find an event in the highest bin). For m^^ — m^^ = 54 GeV, we can combine results 

with the m^^ measurement from WiWi events to deduce that m-^^ ^ 109 ± 3.5 GeV (where 
we have neglected the mismeasurement uncertainty on m{i~^i~)). 

Finally, we turn our attention to top-squarks. In this case, ii — >• bWi, so that stop pair 
production is signalled by events with two 6-jets together with additional jets or leptons from 
the decay products of the charginos and I^j,. Since the titi cross section varies hardly at all 
with PL{e~), we run with right polarized beams (Piie^) = —0.9) to minimize backgrounds 
from WW production. We search for events with > 5 jets, with two tagged as 6's and no 
isolated leptons. We exclude hadronically decaying rs by vetoing jets with one or three 
charged particles as discussed in Sec. 1. This veto capability is crucial to reduce large 
backgrounds from top quark production, where one of the tops decays hadronically and the 
other decays via t hrvr- We also require nji > 140 GeV. For a 50 fb~^ sample of data, 
we are left with a SUSY signal of 286 events, compared with SM background of 36 events. 
At this point, we can plot the energy distribution of the 6-jets, which should have endpoints 
determined by m^-^ and . Again, we fit a function depending on these masses to a large 
sample of generated top-squark pair events, and then obtain contours of Ax^ for 50 fb~^ of 
data. The results are plotted in Fig. 22 a while the resulting Ef, distribution of signal plus 
background is shown in Fig. 22b, along with the best fit histogram. We see that a stop 
mass measurement of = 182 ± 12 GeV is obtained, while = 113.6 ± 8 GeV. This 

measurement of w^fy^ is independent of the measurement from WiWi production described 
above and can provide checks of the inferred sources of the signals. Since fn^^ is determined 

to greater precision via the WiWi channel, the precision of the stop mass measurement can 
be improved by combining the two chargino mass measurements. Since we had used this 
parameter point also to obtain the theoretical fit to the Eb distribution, we have repeated 
this exercise for a somewhat different input for the ti mass in Fig. 21c and Fig. 21d. We see 
that a similar precision is obtained. 



IV. CONCLUDING REMARKS 



Electron-positron collisions provide not only a clean facility for the discovery of super- 
symmetric particles but also provide a unique locale for detailed determination of their 
properties PJTHl^. The discovery of charged sparticles is possible almost all the way up 
to the kinematic limit for their production unless the mass difference between the parent 
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particle and the LSP becomes very small, so that the visible decay products become very 
soft. Their main disadvantages, relative to hadron colliders, are the lower centre of mass 
energy and generally smaller cross sections so that considerable luminosity is needed for 
physics. These are balanced by the clean experimental environment, simplicity of the initial 
state and, at future linear colliders, the availability of longitudinally polarized beams. The 
differences make for complementary capabilities of hadron and e~^e~ colliders. 

In this paper, we have extended the pioneering work of Tsukamoto et. al. and 



analysed what data from experiments at NLC500 might look like if supersymmetry manifests 
itself via the minimal SUGRA framework, and sparticles are kinematically accessible. In 
Sec. II, we have mapped out the discovery reach of NLC500 in the mo— mi/2 plane for several 
sets of other parameters. We have incorporated all the cascade decays of the sparticles into 
our analyses. For the most part, the portion of the SUGRA parameter space that can be 
probed at the NLC is determined by where and Wi signals are observable, though there 
is a small additional region that might be probed via signals from e"'"e~ Z1Z2 production. 
While the availability of longitudinally polarized beams is extremely useful for reducing 
SM backgrounds (mainly from WW production) as well as for isolating various sparticle 
reactions, we found that the beam polarization does not significantly increase the SUSY 
reach if we use the Sex level as our criterion for observability. 

The combined reach of NLC500 (assuming an integrated luminosity of 20 fb~^), in the 
rriQ — mi/2 plane for Aq = 0, tan (3 = 2 and 10 and both signs of /i is shown by the lower solid 
curve in Fig. 23 a- d. These reach contours are not particularly sensitive to the polarizability of 
the electron beam. The dashed and dashed-dotted curves, respectively, denote the kinematic 
limits for producing cr and Wi at Linear Colliders for three different choices of ^/s. The 
small region where the solid curve extends beyond the "kinematic boundaries" is the region 
that should be explored via the Z1Z2 channel. Finally, the upper solid curve denotes the 
reach of the LHC via the J^rp and multilepton channels (the single lepton channel yields 
the greatest reach over the whole space illustrated in Fig. 23) as obtained in Ref. 0. We 
immediately see that as far as, and only as far as, the supersymmetry reach is concerned, 
the LHC reach would be comparable to that of a linear collider operating at a/s ~ 1500 GeV. 
Of course, at this energy SM backgrounds from 2 — 3 and 2 —* 4 (vector boson scattering) 
processes would be important and need to be analysed before drawing definite conclusions. 
We should mention though that NLC500 is guaranteed to find the lightest neutral Higgs 
boson of the minimal model, and so can exclude this framework if no signal is found. In 
contrast, the discovery of the Higgs boson at the LHC is very difficult, and will certainly 
require machine and detector performance at their design levels 

It might also be worth noting that sparticles cannot be much beyond the weak scale if 
SUSY is the new physics that stabilizes the symmetry breaking sector of the SM. Several 
authors |]34| , |35| have attempted to quantify this and obtained upper bounds on sparticle 
masses; e.g. Anderson and Castano have argued that the most favoured region from 
this point of view is where mi/2 < 200 GeV, mo < 200 — 300 GeV. Interestingly, the lightest 
neutralino is an acceptable mixed dark matter candidate if SUGRA parameters are in this 
range ||30|. While this region can partially be probed even at upgrades of the Tevatron [Q, 
and certainly at NLC500, it is worth keeping in mind that fine-tuning considerations are 
qualitative, while the cosmological constraints can be simply evaded, for instance, by allow- 
ing a small amount of i?-parity violation which could have no impact for collider searches. 
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Thus, the larger reach of the LHC or energy upgrades of the NLC provide a safety margin 
and may prove essential for a conclusive exploration of SUSY. Nonetheless, the capability 
of NLC500 for discovering the Higgs boson(s), in itself, appears to us sufficient motivation 
for its construction. 

But the Higgs sector aside, there are many other measurements that are possible at 
Linear Colliders that would be very difficult or impossible at the LHC. These include the 
precision measurement of sparticle masses, spins and couplings which can lead to incisive 
tests of the underlying framework |]14|. In Sec. HI, we have performed four case studies 



to assess the prospects for measuring various sparticle masses. Our study extends the 
earlier analysis in that we allow for all accessible sparticle reactions, and attempt to devise 
strategies to measure slepton masses even when these do not decay directly to the LSP. 
Also, we demonstrate, for the first time, that at least for favourable values of parameters, 
experiments at NLC500, with an integrated luminosity of 20-50 fb~^, should be able to 
obtain masses of sneutrinos, scalar top quarks and even the second lightest neutralino with 
a precision of a few percent. In this analysis, the availability of a polarized electron beam 
with 95% longitudinal polarization has been assumed. We have also identified strategies 
(that again make critical use of the polarization capability of the NLC) to isolate various 
SUSY reactions from one another. This separation facilitates the mass measurements, and 
also allows for the measurement of SUSY cross sections which would be a first step in the 
determination of sparticle couplings. While some mass measurements ||36| are indeed possible 
at the LHC, the systematic precision spectroscopy of sparticles appears to be possible only 
at Linear Colliders. 

As emphasized in Ref. |T^, the measurement of sparticle masses will allow stringent 
tests of the assumptions underlying the SUGRA framework, and so provide a window to 
the symmetries of physics at ultra-high scales. It could be that these assumptions will 
ultimately prove to be incorrect. For instance, other models where SUSY is broken at 
a relatively low scale ~ 10 — 100 TeV have been proposed. These models can have 
significantly different mass patterns and can lead to very different phenomenology |^ from 
what we have considered. Sparticle spectroscopy will provide guidance about the mechanism 
of supersymmetry breaking. On a more pragmatic note, information about chargino and 
neutralino masses and couplings obtained from experiments at NLC500 may prove very 
useful in analysing the complicated cascade decay chains that should be present in the LHC 
data sample. Indeed if it appears that NLC500 is due to become operational significantly 
after the LHC, we would advocate archiving the LHC data in a form suitable for subsequent 
reanalysis in light of new knowledge gained from the NLC. 

In summary, we have affirmed that if sparticles are kinematically accessible at NLC500, 
it will not only be possible to detect the signals, but it will also be possible to measure 
their masses with a precision of 1-5%, even if they do not directly decay into the LSP. 
While electron beam polarization is not essential for SUSY discovery, it is a crucial tool 



14 1 for these precision measurements. The complementary capabilities of NLC500 and 



the LHC cannot be overemphasized PD]. While the LHC essentially probes the complete 



parameter space of the minimal SUGRA framework (mainly via signals from the cascade 
decays of gluinos and squarks), the reach of NLC500 (which indeed probes much of the 
theoretically favoured region) is somewhat smaller. Higgs boson searches, on the other 
hand, are much simpler at the NLC. Also, precision measurements of sparticle masses and 
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couplings, which are only possible at the NLC, would be very helpful in disentangling LHC 
signals from heavy squarks and gluinos (which may not even be accessible at NLC500). 
Working in tandem, experiments at these facilities may allow us to uncover the mechanism 
of electroweak symmetry breaking, and perhaps also provide clues about the symmetries of 
physics at ultra-high energy scales. 
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Appendix: Production Cross-sections with Polarized Beams 



We present in this appendix various SM and SUSY cross sections retaining information 
on the polarization of the incoming beams. The notation used is that of Ref. |T3. First, we 
present lowest order SM cross sections for R or L polarized incoming electrons and positrons. 
For SM fermion and gauge boson pair-production, we have: 



da 

— [eBeL 

dZ L R 



ff) 



An E fL^^' 



where z = cos 9, {6 is the angle between incoming and outgoing fermions) p and E are the 



momentum magnitude and energy of either final state particle, / 
and: 



A" 



Ur, and q, 



^^r{z) 

•' L 



'^1 + z') 



+ ml(l - z^ 



(«£ ± fie 



M|)2 + M|r| 



2{ae±l3e){s-Ml)qf 



([(aJ + /3j)(E^+pV)± 
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E^il + z^)+ml{ 



l-z' 



± 2pfEpz 



da , 

— — (eflCL 

dz L R 



ZZ) 



+ 



4:71 Sy/s 

AMls 



u{z) ^ t{z) 



t{z) u{z) 



1 1 

+ 



u{z)t{z) ^ \t\z) u^{z)^ 
where s, t{z), and u{z) are the Mandelstam variables. 
da 



eReL W+W-) 
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e p 



dz L R 
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where: 



WWR 



4(Q;e + /3e)^tan2 9w\Dz\ 



Ut{z){p's + 3M^) + 4M^p^s^ 



and 



WWL 



Ut{z) 



+4{ae - /9e)'tan2 Owip^s + 3M^)\Dz\ 



Ut{z) 



+ 



+8(Q;e - (5e)ia.n9wM^ReDz + 16(Q;e - /3e) tan^ ^vkM^p^|Dz 



+2 1 - 2(ae - /5e) tan^v^M^i^eDz 



w 



st{z) t{z) 



where C/r(-z) = u{z)t{z) - M^, and Dz ^ (s - + iMzTz)~^. 

The expressions for lowest order MSSM Higgs Production include: 



da 

dz ^"'l"r 



P 



sin2(a + /?) 



167r-y/s sin 6'vk cos^ 6w 



For production, replace sin^(Q; + /3) with cos^(q; + 



— — (e«eL 
dz L R 
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X 



(. - M|)2 + M|r| 

For -ff/i-Hp production, replace cos^(q; + /3) with sin^(Q; + 

— (eRCL H+H-) = —-^(1 - z^) 
dz^ L R ^ 47rVs 



(1 - z^) 



X 



1 



^ 2sin=^gH.- 1 y 
2cos^M/ sin^vK/ (s - M|)2 + M|rl 



+ - 
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For sfermion pair production {fifi, with f — /i, r, i^^u, 1/^-, u, d, c, s, b and i = L or 
R), we find: 



where 
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8q} ^ 2Al{a,± P,fs - 8{a,± Pe)qfAfXs - M^) 



+ 



where Aj^^ = 2(q;/ ^ For the special case of i-ii-i production, we have At^ = 2{af — 
Pf) cos^ 9t + 2(af + Pf) sin^ 9t; for ^2^2 production, simply switch cos^ 9t with sin^ Of. Also, 



da _ ~ 
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For selectron pair production, we find 



where 
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Similarly, 



where 
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where a; 
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For z/e pair production, we find 
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For neutralino pair production, we find 



da 

— ieReL 
dz L R 



ZiZj) 



k 



87rsy/s 



where 



MeeL = 2\Alf\AlfGt{mz^,mz.,me^,z) 



Ae'mj\\ae±Pe)' 
{s - M|)2 + M|r| 



2/2 2 \2 

s -{m^^-m^J 



22 



M: 



ZeR 



_e(-i)(^.+^.+i)(a^ + /?^)(,-M|) 
2[(.-M|)2 + M|r| J 
X [Re ( Wij 5|* ) (m^^ , m , mg^ , z) 



and 



ZeL 



-e(ae-/3e)(g-M|) 

2[(. - M|)2 + M|r|] 

X [i?e ( M/^j Al^. ) (m^^ , m^^. , mg^ , 2;) 



The functions Gi and Gst are defined in Ref. [18 
For chargino pairs, we have 
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and 



M 



ZvL 



8e^(ae-/5e)cot6'vFsin2 7R s(s - M|) 



sin^ Q 



w 



[s - Miy + M|r| 



X 



(xc - yc)[{E - pzf + m|^J + 22/em|,^ 
^2 _|_ p2 _ 2£'p2; + m? 



For W2W2 production, replace Xc with x^, ?/c with y^, sin7/j with cos7/j and with m^^j).^. 
Finally, 
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The final cross section can be calculated from 
o- = fLie~)fLie^)aLL + /L(e")/fl(e+)o-Lij + fB{e~)fL{e^)aRL + fR{e~)fR{e'^)aRR, 

where and /r are defined in Sec. I, and aij = L,R) refers to the cross section from 
e~e^ annihilation. The above formulae have been incorporated into the event generator 
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FIGURES 



FIG. 1. Regions of the uiq vs. mi/2 plane where various sparticle pair production reactions are 

kincmatically accessible. For ZiZj pairs, instead we plot the 10 fb cross section contours. We take 
Aq = 0, tan [3 = 2 and mt = 180 GeV. In a), we take fi < 0, while in b) we take fi > 0. The regions 
denoted by TH are excluded by theoretical constraints, while the region labelled EX is excluded 
by experimental constraints. 

FIG. 2. Same as Fig. 1, except for tan/3 = 10. 

FIG. 3. Various lowest order SM cross sections (in fb) for NLC at ^/s = 500 GeV, as a function 
of e" polarization parameter Pf,(e~). 

FIG. 4. Regions of the mo vs. mi/2 plane where selectrons, charginos and neutralinos are 
accessible at the 5cr level above SM backgrounds, after various cuts discussed in the text. The 

Z1Z2 contour denotes the boundary of the added region where > 10 bb + events (10 fb cross 
section) should be obtained with Pi[e^) = 0.9 above negligible background. We have taken 
^/s = 500 GeV and have assumed 20 fb~^ of integrated luminosity. The SUGRA parameters taken 
are the same as for Fig. 1. In frame a), we show the reach for both polarized and unpolarized 
beams, and in addition, the location of case study points 1, 2 and 3. 

FIG. 5. Same as Fig. 4, except tan /3 = 10. 

FIG. 6. Various lowest order SUSY cross sections (in fb) for NLC at -y/s = 500 GeV, as a 
function of e" polarization parameter PL{e~), for case study point #1. We also show on the left 
the masses of only the accessible superpartners and Higgs bosons. 

FIG. 7. a) Histograms of events after cuts versus missing mass for chargino pair events and 
also for SM and SUSY backgrounds. In b), we plot the Ejj distribution after requiring > 240 
GeV; the arrows denote the kinematic endpoints which are related to the Wi and Zi masses. 

FIG. 8. In a), we plot the distribution after cuts to isolate the Z1Z2 bb + ^ signal from 
SM and SUSY backgrounds. The arrows denote the expected kinematic endpoints of the signal 
distribution. After a cut oirr^> 300 GeV, we plot rrijj for the two detected 6-jets; the distribution is 

related to the light Higgs boson mass, denoted by the arrow. In c), we plot the Ejj distribution 
from the two 6-jets; the kinematic endpoints, which depend on m-^^, m-^^ and on mn^, are denoted 
by the tips of the arrows. 
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FIG. 9. In a), comparison of data to a theory fit depending on and yields the minimum 
value shown, and also the Ax^ = 2.3 and 4.6 contours. In b), a 50 fb^^ 'data' sample is shown 

as well as the best fit theory distribution, for the distribution for Z1Z2 search using case study 

point 1. 

FIG. 10. Same as Fig. 6, except for case study point #2. 

FIG. 11. For 95% right polarized e~ beams, we show in a) the En distribution after cuts for 
case study point #2. The endpoints, denoted by the tips of the arrows, depend on m^^ and rn^^. 
In b), we show the Eg+ distribution which has two main components from eaeR production and 
from crcl production. The corresponding endpoints are shown as well. In c), the E^,- distribution 
has a different structure than that shown in frame b). Finally, after requiring E^+ > 75 GeV in d), 
the Eg- distribution from eRCL is cleanly isolated as discussed in Sec. IIIB of the text. 

FIG. 12. The E^^^""" distribution from Z1Z2 production (after requiring 
along with SM and SUSY backgrounds. The upper endpoint of the solid histogram yields infor- 
mation on m^^ , given knowledge of m^^ and mp,^ . 

FIG. 13. Same as Fig. 6, except for case study point #3. 

FIG. 14. Distribution in E^- from eRSji signal and SM and SUSY backgrounds, for case study 
point #3, using a right polarized e~ beam. The endpoints yield information on mg^ and th-^^. 

FIG. 15. The distribution of Eg in ee/x+jets+^ events from z/gt'e production yields information 
on mp^ and . After comparing 'data' to a theoretical fit, we plot the location of the minimum 

value, and contours of A^^ = 2.3 and 4.6 in frame a). In b), we show the data compared to the 
best theory fit. In c), we again plot A^^ contours for a nearby point in parameter space. Similar 
precision to frame a) is attained. In d), we show the Eg distribution for data and best fit for the 
point taken in frame c). 

FIG. 16. The E^ distribution from the two hardest leptons in slcl — 6£ events leads to a mass 
measurement of mg^ and m-^^. In a), we plot the minimum and also contours of Ax^ = 2.3 and 
4.6. In b), we show a sample of 'data' and also the best fit obtained from frame a). 

FIG. 17. Distribution in dimuon invariant mass showing the Z1Z2 jJ-'^ jJ-' + ^ contribution 
and SM and SUSY backgrounds. The upper endpoint of the signal distribution is bounded by 
m-x — m-x , and leads to a measurement of mx . 

FIG. 18. Same as Fig. 6, except for case study point #4. 
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FIG. 19. In frame a), we show the Ejj distribution from WiWi 1 + jj + If events along with 
SM and SUSY backgrounds. No distinctive endpoints are evident, due to the three-body chargino 
decay kinematics. In 6), we show a scatter plot of Ejj vs. rrijj. The signal is kinematically 
constrained to lie below the solid contour. 

FIG. 20. By requiring rrijj to lie within 4 narrow mass bins, we force 2-body kinematics onto 
the signal distribution. Then, for each bin, the endpoints in the Ejj distribution depend on 
and m~ . In a), we plot the minimum value, along with contours of Ax^, after performing a 
common fit to the four rrijj bins mentioned in the text. In frame 6), we show the 'data' sample 
from the rrijj = 30 GeV bin, compared to the best fit theory distribution. 

FIG. 21. Z2Z2 production leads to a clean sample of e+e^^+//^ + ^ events. A plot of m^+i- 
yields a distribution bounded by m~ — m~ , which gives an estimate of m~ . 

FIG. 22. titi production leads to events containing > 5 jets-|-^, where two of the jets are 
tagged as 6-jets. A plot of Ei, yields information on m^^ and . In frame a), we plot minimum 

and contours of Ax'^ from a comparison of theory to data. In frame 6), we show the 'data' 
sample, and the best fit theory distribution. In c), a similar analysis is made for a nearby point in 
parameter space, and similar precision is obtained. Frame d) shows the associated best fit to the 
'data' distribution. 

FIG. 23. The mo vs. plane is illustrated for Aq = 0, and for tan/? = 2 and 10, and fi < 
and /U > 0. The regions denoted by TH and EX are excluded by theortical and experimental 
constraints, respectively. The dashed lines show contours of mg^ while the dot-dashed lines are 
contours of , and roughly denote the kinematic reach of NLC500, NLCIOOO and NLC1500 
via selectron and chargino pair production channels. The lower solid contour shows the reach of 
NLC500 in parameter space as obtained using our simulation. The upper solid contour shows the 
reach of LHC with 10 fb~^ of integrated luminosity, in the l£+iets+I^rp channel |^. 
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